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We numerically study classical time evolutions of Kaluza-Klein bubble space-time which has 
negative energy after a decay of vacuum. As the zero energy Witten's bubble space-time, where 
the bubble expands infinitely, the subsequent evolutions of Brill and Horowitz's momentarily static 
initial data show that the bubble will expand in terms of the area. At first glance, this result may 
support Corley and Jacobson's conjecture that the bubble will expand forever as well as the Witten's 
bubble. The irregular signatures, however, can be seen in the behavior of the lapse function in the 
maximal slicing gauge and the divergence of the Kretchman invariant. Since there is no appearance 
of the apparent horizon, we suspect an appearance of a naked singularity as the final fate of this 
space-time. 

PACS numbers:04.50.+h, 04.20.Dw, 04.20.Gz 



I. INTRODUCTION 

It is likely that superstring or M-theory governs the physics of gravity or space-time in higher energy stages [Q. 
Such theories are naturally formulated in the higher dimensions than four. We expect a plausible scenario that such 
a higher dimensional space-time somehow evolves to the stable four dimensional space-time according to the history 
of the Universe The so-called brane world scenario (|] is the most actively being investigated along to this line. 

This scenario is motivated by Horava and Witten's theory || which shows that an eleven dimensional supergravity 
theory on the orbihold R 10 x S 1 /Z^ is related to the ten-dimensional E$ x Eg heterotic string theory. Therein the 
matters are confined to the ten-dimensional space-time (three-brane) and gravitons are propagating in the full eleven 
dimensions. The brane world space-time should be stable. 

Although the brane world scenario may be plausible at the reduction from eleven to ten dimensions, the space-time 
will be still compactified to four dimensions in the normal Calabi-Yau's way. Regarding to these full scenario of the 
compactification, the stability of the space-time becomes the important issue to be investigated. The positive energy 
theorem guarantees the stability of the four dimensional asymptotically flat space-time in the framework of general 
relativity (m. Surprisingly, the existence of the extra dimensions can drastically change the situation. Witten showed 
that the five dimensional Minkowski space-time decays into the so-called Kaluza-Klein (KK) bubble space-time unless 
we assume the existence of the elementary fermion related to supersymmetry , of which existence we can not 
expect generally. This also may indicate that the 'bubble' appears somewhere at the bulk or on the brane in the 
brane world scenario and disturbed the three-brane where we are living. 

The metric of the KK bubble space-time given by Witten is written as 

ds 2 = -r 2 dt 2 + (l - d X 2 + (l - -|) ~V 2 + r 2 cosh 2 tdil 2 , (1.1) 

where the x-direction will be compactified and r > ro, and dQ 2 = d0 2 + sin 2 Qdip 2 . In general case, the metric 
has a conical singularity at ro. However, if we carefully take a periodicity along the x-direction, the metric can be 
regularized. More precisely to see this, we write the metric near r = ?*o as 
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dsj ~ -r 2 dt 2 + 2(r ro) d X 2 + , r ° , dr 2 + r 2 cosh 2 t dtt 2 
r A 2(r-r ) 



r 2 cft 2 + 2r [i? 2 d( — V + di? 2 ] + r 2 cosh 2 tdfl 2 , (1.2) 
L V r / J 



where R = y/r — ro- Then we realize that the period should be set to be \ p — 2-ktq. As one can see, the 'boundary 
of bubble' located at r = r$ expands rapidly like coshi and the space-time does not have naked singularities. Here, 
we remind you that the total energy is zero. We are imaging the boundary of the space-time, r = tq, as the surface 
of the 'bubble'. 

Interestingly Brill and Pfister || gave an initial data which has the negative total energy related to the size of 
the compactified dimension. The space-time with the negative energy may be favor in the aspect of energetics. One 
year later, Brill and Horowitz JlO| gave an initial data in a simple way. (We will briefly review their construction in 
Sec. |fi|.) Contrasted to the 'Witten bubble', their solution has arbitrary negative energy regardless of the size of the 
compactified space. This is too far from our intuition that the negative energy is proportional to the Casimir energy 
due to the boundary effect of the compactified space. Therefore it is difficult to imagine the classical evolution after 
the vacuum decay. 

Corley and Jacobson JllJ] discussed the subsequent evolution of Brill-Horowitz's initial data. They found that the 
positive acceleration of the bubble's surface area for the negative mass bubbles, and they conjectured that KK bubble 
with negative energy cannot collapse. However, their study is not sufficient to conclude the final fate of the bubble, 
as they already mentioned, because they considered only the initial behavior of the time-symmetric data and they 
did not do any dynamical studies. 

In this paper, we report our numerical analysis on this final fate problem of KK bubble, especially of the negative 
mass bubble. We start our numerical simulation from the Brill-Horowitz's initial data, and evolve the space-time 
using the standard Arnowitt-Deser-Misner formulation (but 4+1 dimensional decomposition). We will show that the 
space-time initially behaves as Corley- Jacobson's analysis, and expands forever, although the acceleration will be 
negative. Despite of the expanding, we will observe the irregular behavior of the curvature invariant. 

This paper is organized as follows. In Sec. O, we give a brief review of Brill- Horowitz's construction of their initial 



data. In Sec. Ill, we describe numerical method and equations. The results of our simulations are shown in in Sec. 
IV. Finally we summarize our results in Sec. M. 



II. BRILL-HOROWITZ'S INITIAL DATA 

In this section we briefly review Brill and Horowitz's argument [jUj . Let us consider an initial slice with K%j = 
in five dimensional vacuum space-times, where is the extrinsic curvature of a four dimensional spacelike hyper- 
surface. In this slice the Hamiltonian constraint equation becomes ^ R = 0, where R is the four dimensional Ricci 
scalar. Here one can easily see that the Euclidean Reissner-Nordstrom metric with imaginary 'charge' iq satisfies 
the Hamiltonian constraint equation, because the 'energy-momentum' tensor of the four dimensional Maxwell field is 
traceless. The metric of the hypersurface is given by 



(4) .g = U{r)d X 2 



dr 2 

W) 



l dVL 2 



(2.1) 



where U(r) = 1 — 2m/r — q 2 /r 2 and r > r + := m + y/ rn 2 + q 2 . In the same way as the previous Witten's example, 
the metric is approximately written as 



U'(r+) 



R l d 



U'(r+) X \ 2 



dR 2 



r 2 + dfl 2 , 



(2.2) 



near r = r + , where R = y/r — r + . To avoid a conical singularity at r = r + , we assume the period Xp — 47r /U'(r + ) = 
2irr 2 _/{rj r — m) along the ^-direction. 

The total energy is evaluated as E = m/2. The m is arbitrary parameter and q determines the size of the 
compactified space. So the total energy can be arbitrary negative. 
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III. FIELD EQUATIONS AND OUR NUMERICAL METHOD 



To know the final fate of the KK bubbles with the negative energy, we study the subsequent time evolution for a 
long time numerically. We apply 4+1 decomposition of space-time along to the Arnowitt-Deser-Misner formulation 
for the actual time integrations. We describe equations and basic numerical techniques in this section. 

The metric of the full space-time is assumed to be 



d s \ = -N(r, tfdt 2 + e 2a ^U(r)d x 2 + e 2b ^U{r)- 1 dr' 2 + rV^dfi 2 , 



(3.1) 



where U(r) = l — 2m/r — q 2 /r 2 , N is the lapse function, and the metric components a, b and c are now time dependent. 
The evolution equations of the four-metric 7^ and the extrinsic curvature Kij become [] 



K) = N{^R) + KKf] - WrrWDjN, 
jij = -2NK ijt 



(3.2) 
(3.3) 



dot denotes the time derivative, and and denote four dimensional Ricci curvature and the covariant 



derivative, respectively. For the reader's convenience, we write down several terms in (3.2) for the metric fl3.1|) as 



2b 
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r z J 



(3.4) 
(3.5) 
(3.6) 



and 



^D X ^D X N = e 



- 2b 'a' 



^D r ^D r N = e- 2b 



2 U 

(N" -b'N'^U +^N'U' 



^D e ^D e N = e- 2b (c' + ^)N'U, 

where a dash denotes the derivative on r. 

We start our simulation from the initial data of Brill-Horowitz's momentarily static solution, such as 



(3.7) 
(3.8) 
(3.9) 



a(r,0) = b(r,0) = c(r,0) = 0, 
K*(r,0) = K;(r ) 0)=Kt(r,0) = 0. 



(3.10) 
(3.11) 



The numerical region is taken as r + < r < r e , where r + := m + y/ m 2 + q 2 is the location of the bubble at the 
initial data and r e is the numerical outer boundary. We stress from the construction that the Kaluza-Klein bubble 
space-time is restricted in r + < r < 00 jn|. We apply the Robin boundary condition at r — r e such as all the 
components fall off as they form an asymptotically flat spacetime. At the inner boundary r — r + , we use the fact 
that both a and b evolve synchronously as we describe in the Appendix, and use both the evolution equation for trA", 



K A A , . A" - ^D l( ^D z N, 



(3.12) 



where we used the Hamiltonian constraint equation, and the momentum constraint equation, 



MDjK* - {A) D r K = 0, 



(3.13) 



1 Here, for simplicity, we tacitly supposed the boundary condition so that the location of the bubble is 'fixed' under the 
variation of the action. As a result we obtain the 5-dimensional vacuum Einstein equation and can show the consistent result 
given in Appendix A. Since the Cauchy development of the initial data cannot cover all region outside the bubble, one may be 
able to consider another boundary conditions, which might be artificial. 
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so as the system evolves properly. 

In order to specify the lapse function, N, we apply both the geod esic s licing condition, N = 1, and the maximal 



slicing condition, K = 0, which equation becomes [directly from Eq. ( |3.12| )] 



(4) AiV = NK i3 K l: >. (3.14) 

This elliptic equation is solved using the incomplete Cholsky conjugate gradie nt me thod. The outer boundary for N 
is set again as asymptotically flat, and the inner boundary at r + for solving ( |3.14 ) we lineally extrapolate 4-metric 
components. 

We apply Brailovskaya integration scheme (a second order predictor-corrector method) jL4| for the time evolution. 
The accuracy of the calculation is checked by monitoring the violation of the Hamiltonian constraint equation. The 
numerical code passed convergence tests, and the results shown in this paper are all obtained with acceptable accuracy. 



IV. RESULTS 



A. Acceleration of the bubble surface 

We first check whether our code reveals the initial behavior discussed by Corley and Jacobson [jll|. We calculate 
the area of the bubble, 

A(t) = 4ng 0e (r + ,t), (4.1) 

together with its time derivative A, and its acceleration, 

A = Airggg = inr 2 [-2N k e e - 2NK e g + AN 2 {K e g ) 2 ]e 2c . (4.2) 

We first show this acceleration in Fig.[y(a)(b), since this was the quantity discussed by Corley and Jacobson. The 
Figs. 0(a) and (b) are of the geodesic slicing condition and of the maximal slicing condition, respectively. We fix the 
charge q and varied m from negative to positive values. Except for the transition at m = 0, as one can see later, our 
result is not qualitatively sensitive under changes of m/q. Under both slicing conditions, we see that the negative 
mass bubble start expanding (positive A) initially, yet will soon be in de-accelerating phase (negative A), while the 
positive mass bubble keep accelerating all the way in Figs. 0(a), and in the region in Figs. 0(b). More precisely, for 
the positive mass cases in Figs. 0(b), we observe from the numerical results that the acceleration will reach and stay 
at a positive value in the final stage, even if it goes negative for a short time, which is happen to quite small positive 
mass cases. 

Such an initial behavior (for both positive and negative mass bubbles) does agree with Corley and Jacobson's 
analysis (we remark that their analysis was under the geodesic slicing condition). However, the turning behavior into 
de-accelarating phase could not find in their analysis. The de-accelarating does not mean collapsing feature directly. 
Actually, upto we stop our time evolution, the numerical data of the area, (4.1), monotonically increases [Figs. 0(c)], 



while its velocity goes down for negative mass bubbles [Figs. 0(d)]. However, from this facts, we can not say that 
negative mass bubbles will expand forever, because we can see the blow-ups of the Riemann invariant and collapsing 
lapse behavior as we show next. (We had to stop time evolution for negative mass bubble case when we face the 
blow-ups of the Riemann invariant.) 



B. Collapse of lapse 

Since we found that the time integration using the maximal slicing condition survives long term time evolution 
than that of the geodesic slicing condition, we will show only the results of the maximal slicing condition hereafter. 

The maximal slicing condition is known as a robust gauge condition for singularity avoidance (or, exactly speaking, 
avoiding the vanishing of the volume elements of the associated Eulerian observers) jn|. This is because the lapse will 
go quite small value in the strong gravitational field. Contrary, we may guess whether the space-time will collapse or 
not by monitoring the lapse function. 

We plot the lapse function, N, in Fig.||. Fig.g(a) is the lapse function at the bubble surface, r = r+, versus time. 
We see the lapse evolves small value for the case of negative mass bubble space-time. The lines end at the time when 
the violation of the constraint equation begin growing. From above standard behavior of the maximal sliced lapse 
functions, we may say that the negative mass bubble space-time is 'collapsing' in some senses. Fig. 0(b) is snapshots 
of iV at several time for the case of negative mass bubble space-time. 
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C. Riemann invariant 



In order to confirm our guess of the 'collapsing' behavior of the negative mass bubble, we calculated the Krctchman 
invariant (Riemann invariant) RijkiR 1 ^ 1 of both 4 and 5-dimensional Riemann tensor. We see both blow up in the 
cases of negative mass bubbles. We plotted a typical behavior of the invariant as a function of r and t, in Fig.[| These 
lines suggest that the possibility of the formation of singularity in the final phase of evolution. 

D. Apparent horizons 

In order to confirm whether a black hole is formed or not in such a case, we check the appearance of the apparent 
horizon. The location of the apparent horizon is given by the position that the expansion rate, of the outgoing null 
geodesic congruence turns into the negative. The appearance of the apparent horizon indicates the existence of the 
event horizon. 

The definition of the apparent horizon might not be unique in our five-dimensional space-time because it depends 
on the dimension of the space-time which the null geodesic congruence runs^. If the null congruence propagating in 
full five dimensional space-time, the expansion rate, is given by 

= («' + 2c' + \ + ^j)^Je- b - (K* + K$ + K$) (4.3) 

where s a — (0, l/^/g rr , 0, 0) is a outer pointing vector in our spatial four metric. On the other hand, if the null is 
confined to non-compactified four dimensions, the expansion rate, ^9+, is given similarly by 

(3)6» + = ( 3 ) V a s Q -Wjf + s a s b WKab 

= 2(c / + -)VUe- b -(K e g +K^. (4.4) 
r v 

We analyzed both ( 3 )# + and ( 4 )# + in our evolving space-time. Surprisingly, in all cases (positive and negative mass 
bubbles), both expansions remain positive definite everywhere as we show an example in Fig^. These suggest us 
no-appearance of apparent horizons. 

V. DISCUSSION 

We numerically studied the dynamical evolution of the Brill and Horowitz's initial data which can have the negative 
energy. As the zero energy Witten's bubble space-time, we show that the 'bubbles' with negative energy will expand 
by mean of area upto the time we stop the simulations. At first glance this result supports Corley and Jacobson's 
conjecture. However, from the facts that the curvature invariant blows up, and no appearance of the apparent horizon, 
we suspect that a formation of a naked singularity as the final fate of Kaluza-Klein negative energy bubble 0. Hence, 
we may have to consider seriously the decay problem from the Kaluza-Klein vacuum to the Witten-type 'bubble' 
space-time. Possible resolution to this may be given by assuming the supersymmetry which may forbid the decay 
or by constructing quantum gravity theory which may smooth out singularities as normally been expected. 

Although the negative mass bubbles are expanding, we obtained the result that the bubble spacetime terminates at 
the singularity. At first glance, they are incompatible, because the naive picture, which the expanding keeps regularity. 
However, the picture may be based on the Raychaudhri-type equation and the equation does not hold in the present 
case. Moreover, the area cannot properly describes whether the system will collapse or not. Properly speaking, we 
need the proper radius from the center which is absent in the present case. 



2 In usual Kaluza-Klein picture, it is natural that the null geodesic congruence runs in four dimensional part. 

3 An anonymous referee of this article pointed out the similarity of the positive and negative mass bubble results. However, 
from the results we obtained, we believe that there are qualitative differences between positive and negative mass cases in their 
dynamical behaviors. (We remark that our simulations are only up to a finite time in order to keep the resolution against the 
expansion of the spacetime.) 
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Finally, we would like to comment on the so-called brane world scenario |2j |4j]. The brane world is motivated by 
the reduction from the M-theory to the E$ x Eg heterotic superstring theory. This reduction drastically changes the 
picture of the reduction Q because 'matters' are confined to the ten dimensions and gravitons are propagating in 
the eleven-dimensions. Here we call the timelike hypersurface, where matters are confined, by 'brane'. A plausible 
history of the compactification are still quite actively discussed recently and we do not reach the consensus at this 
moment. Apart from this compactification scenario, the reduction from ten to four-dimensions follows the well-known 
Kaluza-Klein type or Calabi-Yau compactification. Thus, our present analysis is basically applicable to the space- 
time on the brane, because we supposed the well-known Kaluza-Klein compactification. More precisely, we may say 
that the Witten-type Kaluza-Klein 'bubble' space-time on the 4-brane will be reduced from at least 6 dimensional 
space-time. As was recently reported ||, the effective Einstein equations on the brane are different from the normal 
Einstein equations. Therefore, it might be worth re-asking what is the final fate of Kaluza-Klein bubble if we describe 
the space-time by such a modified Einstein equation when we take a brane world scenario. 
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APPENDIX A: THE DYNAMICAL EQUATIONS AT THE SURFACE OF THE BUBBLE 



In this appendix, we show that a = b at t he location of the bubble, r = r+, during the time evolution. 
The explicit expression of Eqs ( 3^2 ) and ( 3JJ ) are given by 
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(Al) 



(A2) 



(A3) 



At r — r+, we can truncate U, since U = 0. By adding a suffix + for the variables which is evaluated at r = r + , the 
above equations become 



( a+ 
\N+ 



\N + J 



-N+e- 2b + 



~(d+ + 6+ + 2c + )d+ + l e - 2b +U' + N' + 



-N+e- 2b + 
1 

~~N. 



+ 
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(A4) 



(A5) 



Subtracting Eq. (A5) from Eq. (A4), we obtain 

'a + - K ' 
N+ 



-(a+ +b + + 2c+)- 



(A6) 



G 



Therefore we get 



- b 



A f 4 



± _ Ae -(a++b++2c + ) 



(A7) 



where A is a constant. Since the initial conditions, (3.1C) and (3.11), imply a+ = 
at the boundary, r = r + , we can set a + = b + even after the long time integration. 



which implies A = 0. Therefore 



bep-th/9906064 
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[15] 
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Figure Captions 



Fig.l 

Acceleration of the bubble surface, (4.2), versus time. The figures (a) and (b) are of the geodesic slicing condition 
and of the maximal slicing condition, respectively. In both figures, we see that the negative mass bubble will soon be 
in collapse phase, although they start expanding initially. We set q = 1 (hereafter for all figures). 
For the case of evolutions with the maximum slicing condition, we plot (c) area of the bubble surface A, (4.1), versus 
time, and (d) the velocity of the bubble surface, dA/dt. 



Fig.2 

The lapse function, N , are plotted (the solutions of maximal slicing condition). The figure (a) is N at r = r + versus 
time. We see the lapse evolves small value for the case of negative mass bubble space-time. The figure (b) is snapshots 
of N at several times for a case of negative mass bubble (to = —0.4) space-time. 



Fig.3 

(a) Typical snapshots of the Ricmann invariant RijkiK 1 -' of 4-dimensional Riemann curvature for the case of negative 
mass bubble (to = —0.4) are plotted. Only the region near the bubble surface is drawn, (b) The Ricmann invariant 
RijkiR 1 ^ 1 °f 5-dimensional Riemann curvature are plotted as a function of time. We see blow-ups in the cases of 
negative mass bubbles (we cut the display range at 10 5 ). The values are evaluated at a point right from the bubble 
surface (that is, at r + + A?') . 



Fig.4 

A typical sample of the outgoing null expansion ra,te^6 + and*- 4 - 1 ^ are plotted for the case of negative mass (m = —0.4) 
bubbles. 
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FIG. 1. Acceleration of the bubble surface, d 2 A/dt 2 , versus time. The figures (a) and (b) are of the geodesic slicing condition 
and of the maximal slicing condition, respectively. In both figures, we see that the negative mass bubble will soon be in collapse 
phase, although they start expanding initially. We set q = 1 (hereafter for all figures). 
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(c) maximal slicing evolution 




15.0 
time 




FIG. 1. (continued) 

For the case of evolutions with the maximal slicing condition, we plot (c) area of the bubble surface A, ( |4,l| ), versus time, and 
(d) the velocity of the bubble surface, dA/dt. 
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t=0 . 

t = l - 5 

t = 3 . 

t = 4 . 5 

■ - - t=6 - 
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m = + 0.2, 0.4 



0.0 5.0 10.0 15.0 20.0 25.0 30.0 

time 

FIG. 3. (a) Typical snapshots of the Riemann invariant RijuR'^ 1 " 1 of 4-dimensional Riemann curvature for the case of 
negative mass bubble (m = —0.4) are plotted. Only the region near the bubble surface is drawn, (b) The Riemann invariant 
RijkiR 1 ^ 1 of 5-dimensional Riemann curvature are plotted as a function of time. We see blow-ups in the cases of negative 
mass bubbles (we cut the display range at 10 5 ). The values are evaluated at a point right from the bubble surface (that is, at 
r+ + Ar). 
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FIG. 4. A typical sample of the outgoing null expansion rate ^0+ and ^ 4 '0+ are plotted for the case of negative mass 
(m = -0.4) bubbles. 
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